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Abstract 

The aim of the article is to construct the S'-matrix interpretation of the perturbation 
theory for the Wigner functions generating functional at a finite temperature. The tem- 
perature is introduced in the theory by the way typical for the microcanonical description. 
The perturbation theory contains the two-temperature Green functions. The two possi- 
ble boundary conditions are considered. One of them is usual in a field theory vacuum 
boundary condition. The corresponding generating functional can be used in the particle 
physics. Another type of the boundary condition assumes that the system under consider- 
ation is in environment of the free particles background field. This leads to the theory with 
Kubo-Martin-Schwinger boundary conditions ones at the one-tcmpcraturc limit. 
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1 Introduction 



The field-theoretical description of statistical systems at a finite temperature is usually 
based on the formal analogy between imaginary time and inverse temperature /?(/? = 1/T) 
. This approach is fruitful 0] for description of the static properties of a system, but it 
demands a complicated mathematical apparatus for the analytic continuation to the real 
time , if we want to clear up the dynamical aspects. At the same time a great number 
of modern problems is connected with system dynamics. It is important, for instance, for 
description of evolution of the Universe, of quark-gluon plasma behavior, etc. 

The first important quantitative attempt to build the real-time finite-temperature field 
theory discover a problem of the pinch-singularities. The further investigation of the 
theory have allowed to demonstrate the cancellation mechanism of these unphysical singu- 
larities ||. This attained by doubling of the degrees of freedom |], [7|: the Green functions 
of the theory represent 2x2 matrix. It surely makes the theory more complicated, but 
the operator formalism of the thermo-field dynamics §| shows the unavoidable character 
of this complication. 

It is important to note that the traditional real-time finite-temperature field-theoretical 
description || [7]] of statistical systems based on the Kubo-Martin-Schwinger (KMS) || 
H condition for a field: 

= (1.1) 

which, without fail, leads to the equilibrium fluctuation-dissipation conditions JO]] (see 
also pp. 

We shell use the S"-matrix approach which is natural for description of time evolution. 
(The S'-matrix description was used also in El |T^|.) For this purpose the amplitudes 



< {p)m\{q)n >= a n>m (pi,p 2 , -;P m ; qi,q 2 , -An) (1-2) 

of the n- into m-particles transition will be introduced (Sec. 2). The in- and out-states 
must be composed from mass-shell particles fT5 |. Using these amplitudes we will calculate 



the probability which is 

~ \a n ,m\ 2 =< (P)m\(q)n >< {q)n\{p)m > ■ (1-3) 

This will lead to the doubling of degrees of freedom. 

The finite temperature will be introduced (see also [[16]) taking into account that, for 
instance, 

^-^'fiefe (1.4, 

is the differential measure of the initial state. Then we will define the temperature as the 
function of initial energy through the equation of state (see Sec. 2). This introduction of 
temperatures as the Lagrange multiplier is obvious for microcanonical description ||. 

The final-state temperature will be introduced by the same way. So, we will construct 
the two-temperature theory. In this theory with two temperatures (the equation of state 



can be applied at the very end of calculations) it is impossible to use the KMS boundary 
condition. 

One can note that the product of amplitudes ( |1.3|) describes the closed-path motion in 
the functional space, with some "turning-point" fields $(o"oo), where o"oo is the infinitely 
far hypersurface (Sec. 3). The value of $(o"oo) specify the evironment of a system. 

In Sec. 2 the vacuum boundary condition ^(cToo) = 0, familiar for a field theory, will be 
considered. This theory can be applied in the particle physics. The simplest (minimal) 
choice of ^(a^) ^ assumes that the system under consideration is surrounded by black- 
body radiation (Sec. 4). This interpretation restores Niemi-Semenoff's formulation of the 
theory 0. 

One should admit also that this choice of boundary condition is not unique: one can 
consider another organization of the environment of considered system. The S'-matrix 
interpretation will be used since it able to show the way of adoption of the formalism to 
the arbitrary boundary conditions. It should broaden the potentialities of the real-time 
finite-temperature field-theoretical methods. The special interest represent the topological 
effects, but in this paper consideration will be performed in the perturbation framework 
only. 



2 Vacuum boundary condition 

The starting point of our calculations is n- into m-particles transition amplitude a ntTn , 
the derivation of which is well known procedure in the perturbation framework. Let us 



introduce (n + m)-point Green function G n ^ m through a generating functional Zj |L7 



G n , m ((x) n , (y) m ) = Hr +m n jm n KvjZj, (2.i) 



5j{x) 

Zj = [ D$>e lS ^\ (2.3) 



where 

and 
The action 

Sj($) = - V($) + J dxj(x)$(x), (2.4) 

where S($>) is the free part and V($) describes the interactions. At the end one can put 
3 — 0- 

To provide the convergence of the integral ( |2.3| ) over scalar field $ the action Sj(§) 
must contain positive imaginary part. Usually for this purpose Feynman's ie-prescription 
is used. It is better for us to shift infinitesimally the time contour to the upper half plane 



18], 13], i.e. on the following contour C + : t — > t + ie, e > and after all calculations to 



return the time contour on the real axis, e — > +0. 



In eq. ( |2.3| ) the integration is performed over all field configurations with standard 
vacuum boundary condition: 

/ d^xd^d^) = [ da^d^ = 0, (2.5) 

which presides to zero contribution from the surface term. 
Let us introduce now field (f> through the equation: 

and perform the shift $ — > $ + <fi in integral (|2.3| ), conserving boundary condition (|2.5|) . 
Considering cf) as the probe field created by the source: 

(j){x) = J dyG (x-y)j(y), 
{d 2 + m 2 ) x G (x-y) = 6(x-y), (2.7) 

the only connected Green function G c n m will be interesting for us. Therefore, 

n m 

G c n , m ((x) n , (y) m ) = (-t) n+m II fa k ) II Kvk)Z{4>), (2-8) 

fc=i fc=i 

where 

Z{<f>) = J D^e iS{ -^- iV ^ +4>) (2.9) 

is the new generating functional. 

To calculate the nontrivial elements of S-matrix we must put the external particles on 
the mass shell. Formally this procedure means amputation of the external legs of G c nm 
and further multiplication on the free particles wave functions. In result the amplitude of 
n- into m-particles transition a„ im in the momentum representation has the form: 



n 



a n ,m((l)n, (P)m) = H)^™ II <K%) II P (j>k)Z{(l>) . (2.10) 

fc=l fc=l 

Here we introduce the operator 

0(g) = J dxe-^fc), 4>{x) = j—- (2.11) 

and qk and pk are the momentum of in- and out-going particles. 

Supposing that the momentum of particles are insufficient for us the probability of n- 
into m-particles transition is defined by the integral: 



n\m\ 



/n m 
du n (q)du m (p)5 {4) {^2 1k~Yj 
Pk) \0"n,m \ , (2-12) 

k=l k=l 



where 

du n ( q ) = n Mo*) = n n tl k ( v e = + m2 ) 1/2 ' ( 2 - 13 ) 

fc x = i k=i (27r) 3 2e(g fc ) 

is the Lorentz- invariant phase space element. We assume that the energy-momentum 
conservation 5-function was extracted from the amplitude. 

Note that r n m is the divergent quantity. To avoid this problem with trivial divergence 
let us divide the energy-momentum fixing ^-function into two parts: 

<5 (4) (E^-E^) = / d A P5^\P-Y.q k )t (A \P-Y,V k ) (2.14) 
and consider a new quantity: 

i „ n n 

r ( P ) = Y,-n / dco n (q)dw m (p)5^\P-J2Qk)S^\P-^Pk)K m \ 2 - (2-15) 

— 11 'ill ' I — — 

n,m IL -" L - J k=l k=l 

Here we suppose that the number of particles are not fixed. It is not too hard to see that, 
up to phase space volume, 

r = J d 4 Pr(P) (2.16) 

is the imaginary part of amplitude < vac\vac >. Therefore, computing r(P) the standard 
renormalization procedure can be applied and the new divergences will not arise in our 
formalism. 

The Fourier transformation of ^-functions in ( |2.15| ) allows to write r(P) in the form: 



(2vr) 4 (2tt) 



where 



i „ n m 

%,02) = E -Tl / UiMqk)e- ia ^} l[{MPk)e- ta2Pk }\a n , m \ 2 . (2.18) 

— n'm' I 

n,m •> k = i k = i 

The introduction of the "Fourier-transformed" probability -R(ai,a:2) means only that the 
phase-space volume is not fixed exactly, i.e. it is proposed that 4-vector P is fixed with 
some accuracy if «j are fixed. The energy and momentum in our approach are still locally 
conserved quantities since an amplitude a nm is translational invariant. So, we can perform 
the transformation: 

«iE^ = ( a i -^lEft + ^E* ( a i ~~ E + °\P (2-i9) 

since 4-momenta are conserved. The choice of <Ji fixes the reference frame. This degree 



of freedom of the theory was considered in |L9[ ^0 
Inserting (|2.10|) into (|2.18|) we find that 



i?(«i,a 2 ) = exp{i J dxdx'((f) + (x)D + _(x — x', a 2 )0_(x') — 

4>4x)D„ + (x - x', ai )Mx'))}Z^ + )Z*(^), (2.20) 



where D + _ and D_ + are the positive and negative frequency correlation functions corre- 
spondingly: 

D + _(x - x', a) = -i I du{q)e iqix - x '- a) (2.21) 



describes the process of particles creation at the time moment xq and its absorption at 
Xq, Xq ^> Xq, and a is the center of mass (CM) 4-coordinate. Function 

D_ + (x -x',a)=i J duj(q)e- iq(x - x ' +a) (2.22) 

describes the opposite process, xq < x' Q . These functions obey the homogeneous equations: 

(d 2 + m 2 ) x G + ^ = (d 2 + m 2 ) x G^ + = (2.23) 

since the propagation of mass-shell particles is described. 

We suppose that Z(<p) may be computed perturbatively. For this purpose the following 
transformation will be used: 

e -iV(4>) _ e ~i J dxj{x)4>'(x) e i J dxj{x)<j>{x) e -iV{4>') = 

e J dx<f>(x)j>'(x) e -iV (<!>') 



-iV{-i]) e i J dxj(x)<f>{x) ^ (2.24) 



where j was defined in ( p.2|) and <j> in (|2.11|) . At the end of calculations the auxiliary 



variables j, <p' can be taken equal to zero. Using the first equality in ( |2.24| ) we find that 



Z((J)) = e~* f dx Kx)Hx) e -iV($+'t>) e ~{ J dxdx 1 j(x)D ++ (x-x')j(x') ^ (2.25) 

where D ++ is the causal Green function: 

(d 2 + m 2 ) x G ++ (x -y) = 5(x - y) (2.26) 

Inserting ( |2.25| ) into ( |2.2(J| ) after simple manipulations with differential operators, see 
( |2.24| ) we find the expression: 

R(a u a 2 ) = e -m-i'U)+iV(-ij-) x 

x exp{^ J dxdx'{j + {x)D + ^(x - x', a^j^x') - j^{x)D_ + (x - x', a 2 )j+(x') - 

-j + (x)D ++ (x - x') U (x')+j4x)D_4x - x')j4x'))}, (2.27) 

where 

= (D ++ y (2.28) 

is the anticausal Green function. One can consider i?(ai,«2) = R{oti, 012) 31,32) as the 
generating functional for Wigner functions (see also |3], |2~T|]). 

Considering the system with large number of particles we can simplify calculations 
choosing the CM frame P = (Pq = E,0). It is useful also |l6|, || to rotate the contours 



of integration over ao,k- a o,k = —i(3k, Imfik — 0, k = 1, 2. In result, omitting unnecessary- 
constant, we will consider R = R(p 1 ,p 2 )- 

External particles play the double role in the ^-matrix approach: their interactions 
create and annihilate the system under consideration and, on the other hand, they are 
probes through which the measurement of a system is performed. Since /3k are the conju- 
gate to the particles energies quantities we will interpret them as the inverse temperatures 
in the initial and final (Z^) states of interacting fields. But there is the question: are 
constants (3% really the "good" parameters to describe the system. 

The integrals over (3k- 

r f E ) = I ^^(fc+A^e-^A), (2.29) 

where 

F(J3 1 ,p 2 ) = -hiR(p 1 ,p 2 ), (2.30) 

can be computed by the stationary phase method. This assumes that the total energy E 
is a fixed quantity. The solutions of the equations (of state): 

E = m^M (231) 

o/3k 

gives the mostly probable values of (3f. at a given E. Eqs. ( |2.31| ) always have the real 
solutions and, because of energy conservation law, both eqs. (|2.31|) have the same solution 
with the property ||: 

Pk = (3(E), (3>0. (2.32) 

If the fluctuations of (3k are large it is insufficient to know (3(E): the expansion of inte- 
gral (|2.29|) over ((3 — (3k) will lead to the asymptotic series with zero convergence radius 
since F(/3x,/3 2 ) is essentially nonlinear function. In this paper we assume that (3 is the 
"good" parameter, i.e. the fluctuations of (3k are Gaussian. In this case we can interpret 
F{(3i,(32) as the free energy and l/(3k as the temperatures. Such definition of thermody- 
namical parameters is in a spirit of microcanonical description. Note the important role 
of decomposition (2.14) in this interpretation of the (3k- 



The structure of generating functional ( |2.27|) is the same as the generating functional of 



Niemi-Semenoff [KJ have. The difference is only in the definition of Green functions which 
follows from the choice of boundary condition ( |2.5|) . The Green functions Dy, i, j = +, — 
were defined on the time contours C± in the complex time plane (C_ = Cl). This 
definition of the time contours coincide with Keldysh' time contour [^j. The expression 
Q2.271 ) can be written in the compact form if the matrix notations are used. Note also a 
doubling of the degrees of freedom. This doubling is unavoidable since Green functions 
Dij are singular on the light cone. 



3 Closed-path boundary conditions 

The calculation of generating functional R{a,i,a,2) is performed introducing the corre- 
sponding generating functional 

Ro{<j> ± ) = Z{<t> + )Z*{<j>_) = 
J exp{z5($+) - *S($_) - iV{§+ + </>+) + + </>_)}, (3.1) 



see ( |2.20| ). The fields and $ + ,$_ were defined on the time contours C + ,C_. By 

definition, path integral (|3.1|) describes the closed path motion in the space of fields $. 
We want to use this fact and introduce a more general boundary condition which also 
guaranties the cancelation of the surface terms in the perturbation framework. We will 
introduce the equality: 



f da^+d^ = f da^-d^-. (3.2) 



The solution of eq. (|3.2[ ) requires that the fields $ + and $_ (and theirs first derivatives 
d fl &±) coincide on the boundary hypersurface a^: 

$ ±(<Toc) = H(Too), (3.3) 

where, by definition, $((7oo) is the arbitrary, "turning-point", field. 

In absence of the surface terms, the existence of nontrivial field ^(a^) has the influence 
only on the structure of Green functions 

G++ =< T$+$+ >, G+_ =< >, 

G_+ =< >, G__ =< f > (3.4) 

where T is the untitemporal time ordering operator. This Green functions must obey the 
equations: 

(d 2 + m 2 ) x G + 4x -y) = (d 2 + m 2 ) x G„ + {x -y)=0, 
(d 2 + m 2 ) x G ++ (x -y) = (d 2 + m 2 )* x G^(x - y) = 5(x - y), (3.5) 

and the general solution of these equations: 

Ga Dii ~\~ Qiii 

Gij = g ij: i ^ j (3.6) 

contain the undefined terms which are the solutions of homogenous equations: 

(d 2 + m 2 ) x g ij (x -y) = 0, i,j = +,-. (3.7) 

The general solution of these equations (they are distinguished by the choice of the time 
contours C±) 

9ij {x -x')= [ duiqy^-^n^q) (3.8) 



are defined through the functions n%j which are the functionals of "turning-point" field 
^(croo): if $(o"oo) = we must have = and we will come back to the theory of 
previous section. 

Our aim is to define n.^. We can suppose that ~< $(o"oo) • ■ ■ ^(^oo) >• The 
simplest supposition gives: 

riij ~< >~< $Voo) > . (3.9) 

We will find the exact definition of starting from the ^-matrix interpretation of the 
theory. 



4 KMS boundary condition 

In the previous section it was shown that the theory permits the arbitrariness of boundary 
condition: the turning-point field $(000) may be arbitrary since the "closed-path" motion 
in the functional space is described. We will suppose that on the infinitely far hypersurface 
(Too there are only free, mass-shell, particles. Formally it follows from (|3.6|) - (|3.8| ). This 
assumption is natural also in the S'-matrix framework f!5 |. In other respects the choice 
of boundary condition is arbitrary. 

Therefore, our aim is concerned with description of evolution of the system in a back- 
ground field of mass-shell particles. In this paper we will assume that there are not any 
special correlations among background particles. We will take into account only the re- 
strictions connected with energy-momentum conservation laws. Quantitatively this means 
that multiplicity distribution of background particles is Poison-like, i.e. is determined by 
the mean multiplicity only. This is in spirit of definition of in eqs. ( |3.8| ), (|3.9| ). 



Our derivation is the same as in p0[ . Here we restrict ourselves mentioning only the 
main quantitative points. 

Calculating the product a„ im a* m we describe a process of particles creation and further 
their adsorption. In the vacuum case of Sec. 2 the two process were taken into account: 
creation of particles by the a niTTt and theirs adsorption by a* m , and the opposite process of 
particles creation by a* m and theirs adsorption by a n , m . This processes were time ordered. 

This was the reason of the frequency correlation functions D + - and D y appearance. In 

the nonvacuum case, i.e. in presence of the background particles, this time-ordered picture 
is slurring over since the possibility to absorb particles before their creation appears. 

The processes of creation and adsorption are described in vacuum by the product of 



operator exponents of </>+</>_, 4>-4>+- We can derive (see also [2D|) the generalizations of 
( |2.2C| ): presence of the background particles will lead to the following generating func- 
tional: 

R cp = e iN ^R (<P ± ), (4.1) 

where Ro(<fr±) is the generating functional for vacuum case, see ( |3.1| ). The operator 
N(4>*(f)j),i,j = +, — , describes the external particles environment. 

The operator <j)*(q) can be considered as the creation and <fii{q) as the annihilation 
operator, see definition ( p. 10 ). Correspondingly the product 0*(<?)0j(<7) acts as the activity 



operator. So, in the expansion of N{4>*<j)j) we can leave only the first nontrivial term: 

N{fcfa) = I tM9)&(?Hife(?), ( 4 -2) 



since no special correlation among background particles should be expected. If the exter- 
nal (nondynamical) correlations are present then the higher powers of <f)*4>j will appear in 
expansion ( f4.2| ). Following to the interpretation of 4>*4>j we conclude that is the mean 
multiplicity of background particles. 
In ( f4.2p the normalization condition: 

N{0) = (4.3) 

was used and summation over all i,j was assumed. In the vacuum case only the combi- 
nations i 7^ j are present since the time ordering. 

Computing R cp we must conserve the translational invariance of amplitudes. It lead 
to extraction of the energy-momentum conservation (^-functions. Having background par- 
ticles flow it is important to note that only connected contributions must be taken into 
account. Using the Fourier transformation we can find that to each vertex of in- going in 
a> n ,m particle we must adjust the factor e~ iaiq ^ 2 and for each out-going particle we have 
correspondingly e~ ia2q / 2 . 

So, the product e ~ takq l 2 e~ ta:iq l 2 can fog interpreted as the probability factor of the 
one-particle {creation + annihilation) process. The n-particles {creation + annihilation) 
process' probability is the simple product of these factors if there is not the special cor- 
relations among background particles. This interpretation is evident in the CM frame 
«fc = {-i(3k,0). 

After this preliminaries it is not too hard to find that in the CM frame we have: 

V°° tip — a 4 9oh 
n++{qo) = n__(g ) = h+02 . , = 

V°° p- 2 ™ \n 
* ~ / 1 I /^l + /?2 n (a a \ 

7^m =w(|qb| ^~ ) - (4 - 4) 

Computing n^ for i ^ j we must take into account that we have one more particle: 

V°° iip- a P V°° np 2 ion 

n+-(9b) = g(ft) ^ + e(-g ) L -°^ 1+fa = 

6(g )(l + n{qM) + e{-q )n{-q Pi) (4.5) 

and 

n„ + (g ) = e(g )n(go/3 2 ) + e(- 9o )(l + n{-qM). (4.6) 
Using ( |4.4| ), ( f4.5|) and ( (4.6[) , and the definition ( |3.6|) we find the Green functions: 

G M (x - x', 09)) = / -^e^^Gdiq, {[3)) (4.7) 



where 



^j(g,(/5)) = r 2 -Q 2+4e + 

\ q 2 — m 2 —ie / 

+ 2^- m >)(*^f g . fi ™ +( * > ) (4-8) 



and 



a± (/5) = _ e §(M±*)). (4.9) 
The corresponding generating functional has the standard form: 

# cp (j±) = exp{-zF(-zj + ) + iV(-ij-)} x 

xexp{^ / dxdx'ji(x)Gij(x — x',(j3))jj(x')} (4.10) 



where the summation over repeated indexes is assumed. 

Inserting (|4.10|) in the equation of state ( |2.31j ) we can find that Pi = P2 = P(E). If 
(3{E) is a "good" parameter then GV,-(x — x'\ 0) coincide with the Green functions of the 
real-time finite-temperature field theory and the KMS boundary condition: 

G + _(t-t') = G_ + (t-t' -i[3), G-+(t-t') = G+-(t-t' + i/3), (4.11) 

is restored. The eq. ( |4.11| ) can be deduced from (|4.8|) by the direct calculations. 



5 Conclusion 

In our interpretation of the real-time finite-temperature field theory the statistics and 
fields quantum dynamics were unlinked: statistics is fixed by the operator exp{iN((f>*(j)j} 
and a pure field-theoretical dynamics is described by Ro((fi±) = Z(4> + )Z*((f)_), where 
Z((p±) is the vacuum into vacuum transition amplitude in presence of the external (auxil- 
iary) fields < vac\vac >^. We can say that the operator exp{iN(<pi(j)j)} maps the system 
of interacting fields on the state with definite thermodinamical parameters. 

This means the following procedure: we can act by the operator exponent expjiV} 
before calculation of Rq. We find in this case from (4.10) that 



R cp (Pi,p2) = J D$ + £$_ e ^(<M-«(<M e iW±) ; (5.1) 
where 5($) is the total action and 

$ ± ) = _, ln { e ^($ + +<M-;v(<i>-+«M 

e i J dYdy[$i(Y+y/2)G ij (yM)fa(Y-ym e -iV(*++cl> + )+iV(*-+^)y z^) 

is act as the effective source of fields $± and depends from the temperature Such 
representation introduce in quantization procedure the external conditions of the consid- 
ered problem. This can help to calculate the observables with high accuracy. 



The formalism allows to describe an arbitrary system. The special interest presents the 
"local equilibrium" case. The structure of R cp and of Green functions remains the same 
for this case, see (|4.10| ) and ( |4.S| ), but [3k becomes coordinate dependent: [3k = PkiX)- 
Note that Y = (x + x')/2 is the Wigner's coordinate [14], 
will be given in the following paper. 



HI! . The derivation of this result 
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